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THEORETICAL STUDY OF THE COUPLING 
BETWEEN THE SOLAR WIND AND THE EXOSPHERE 

APPENDIX A 

Superthermal Electron Production in the Transition Region* 

R. W. Fredricks, F. L. Scarf and W. Bernstein 
TRW Space Technology Laboratories 

The ion acoustic wave -electron gyrofrequency resonance is examined 
quantitatively and the results used to discuss the production of E =» 30 kev 
electrons in the solar wind-geomagnetic field transition region. The wave 
is represented by a potential distribution <|) = (J) q sin(k • r - cut +\jf) with 
0 < u) < Q. (ion), 0 < k < k (Debye), and \jj a stochastic phase variable 
associated with nonlinear effects. It is estimated the e(j) o may approach 
35$ of the incident wind energy in localized regions, and sample numerical 
solutions to the equations of motion are used to show that the maximum 
electron energy derived from this resonance is of the order of 130 e(J) Q . 

The resonance conditions leading to production of significant (30-80) kev 
electron fluxes are discussed and sane favorable magnetic field configura- 
tions are implicitly described. 


* Transactions, A.G.U., 4_5, No. 4, 630, December 1964. 
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THEORETICAL STUDY OF THE COUPLING 
BETWEEN THE SOLAR WIND AND THE EXOSPHERE 

APPENDIX B 

A Possible Interpretation of the P-11 VLF Measurements 

F. L. Scarf and R. W. Fredricks 
TRW Space Technology Laboratories 

The possibility that the large amplitude VLF electric fields detected 
on the P-11 spacecraft are associated with longitudinal plasma oscillations 
is considered. In a collisionless plasma, the mode that can be generated 
most easily is the ion acoustic mode with E = - V(j), <|) = (j) Q cos (k • "r - cut) 
and 0 < <x> «; Q ( ion) . If T =» =» T . , theory indicates that a very small 
current should induce growing waves, but laboratory experiments show that 
the instability occurs even more readily than predicted by the theoretical 
threshold calculations. These electric field fluctuations allow fast 
diffusion of charged particles across magnetic fields and thus the presence 
of ion acoustic waves could be related to electron precipitation and 
diffusion phenomena observed on Injun 3, Explorer 14, and Air Force 
satellites. The mechanism for auroral injection may also involve the fast 
diffusion process. 


* Transactions, A.G.U., 45, No. 4, 598> December 1964. 
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CONDUCTIVE HEATING OF THE SOLAR WIND, II: 

THE INNER CORONA 

N65-20622 

F. L. Scarf and L. M. Noble 


TRW Space Technology Laboratories, Redondo Beach, Calif. 


Abstract 


2 5 k 2. 


New numerical solutions to the Navier -Stokes equations 
are used to describe the streaming solar atmosphere. It is 
argued that no external heat input is needed beyond about 
2R @ if the coronal base temperature is near 1.5 x 10° °K 
and the plasma has a uniform 10# concentration of helium. 

In the subsonic region, thermal conduction dominates the 
flow and viscous effects are unimportant. The transition 
to supersonic flow occurs near 7R @ > an< i beyond this radius 
viscous forces are extremely significant. It is shown that 
solutions to the inviscid model are singular and that they 
are not related to the limiting solutions for small, but 
finite, viscosity. 


Submitted to the Astrophysical Journal, November 1964. 
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1. INTRODUCTION 

Recent numerical calculations (Noble and Scarf, 1963; hereafter referred 
to as I) tend to confirm the hypothesis that the solar wind is a natural ex- 
tension of the solar corona with thermal conduction beyond a thin shell heat- 
ing region providing the main local energy flux which induces the supersonic 
flow. In this series of papers the validity of the conductive heating model 
is examined in greater detail, and the model is modified and extended. The 
previous discussion was based on study of a particular solution to a dimen- 
sionless form of the Navier-Stokes equations, and viscous dissipation was not 
explicitly included. These calculations have now been supplemented in two 
ways: (a) a broad range of solutions is examined in order to determine the 

best fit to the observed coronal density in the inner region, R^"* r < (8-10)R @ . 
As in (i), no external heat source is considered and viscous terms are neg- 
lected near the base; (b) viscous corrections are inserted and it is found 
that they become significant beyond the sonic transition at r = r so that 
subsequent numerical integration must be performed without neglect of these 
terms. A brief summary of some of these results was recently presented 
(Scarf and Noble, 1964), and in this paper we wish to describe these refine- 
ments in detail. 

It will be shown that the inviscid model represents a singular flow 
pattern, and that the correct flow exhibits extremely large velocity gradients 
so that the applicability of the Navier-Stokes equations becomes questionable 
at large distances ^(15 - 20 )R q J . The final paper in this series will contain 
an investigation of the breakdown of the Navier-Stokes relations, together with 
speculation on the transition to a form of modified free molecular flow. 


-la- 


2. THE CONTINUUM EQUATIONS 


The solar corona is a fully ionized gas containing approximately 10$ 
helium and 90$ hydrogen. This plasma is described by a series of Boltzmann 
equations for electron, proton and alpha -particle distribution functions, 
and the equations are strongly coupled by coulomb collision terms. If all 
species have the same local temperature, T, and the same streaming velocity, 
u, (i.e., if diffusion effects can be neglected), then the lowest. order 
mojnents of the distribution functions satisfy the following general equations 
(Burgers, i960). 


ff + ^ (Nu i ) = ° ' « 

H ^ + N \ 3^ + m ^ (NkT) + m ^ P ik ‘ Bs i = 0 ' < 2 > 

lit <"> + l“i * I ** 3^ + \ P ik € ik + - °’ <3) 


with 



(4) 


Here N is the total number density, (N + N + N ), m is the mean mass per 

c Jr ^ 

particle, Nmg^ is the external force per unit volume, and Q^, P represent 
sums over all species of higher moments of the distribution functions, 

“i (£><•) = I f \i - sto) 2 ■ a(E»i f s te t > a3v - (?) 

s J 

e,p,a 

F ikt' t) - £ m s 

s 

kT^tjS^ 



- 2 - 


( 6 ) 


The coupled Boltzmann equations lead to very complex general relations 

for P ik (Burgers, i960), "but these can be considerably simplified if each 

f (v,r,t) is sufficiently close to a stationary local Maxwell -Boltzmann dis- 
3 ** *** 

tribution. 


,( 0 ) 


r m 

[25%. 


13/2 


exp 


* 

-m (v - u(r))* 
2kT(r) 


■} 


(7) 


and if the gradients of N, T and u are "sufficiently small". Although it 
will be necessary to return to this definition of "sufficiently small" in a 
subsequent discussion, we presently assume that in the lower corona the afore- 
mentioned conditions are satisfied so that the conventional linear approxi- 
mations, Q » -WT, are valid; Eqs. (l) - (3) are then the 

Navier -Stokes equations for motion of a fluid with a mass density 

Nm = Zn m and a specific heat ratio of 5/3. 
s s 

The transport coefficients ~K, are evaluated by examining the colli- 
sion terms in the coupled Boltzmann equations. Chapman (1954) computed X 
and 'ty in this way for a fully ionized hydrogen gas, and he noted that the 
electrons are primarily responsible for thermal conduction while the ions 
dominate the viscous contribution. For a hydrogen plasma, the Chapman 
transport coefficients are 

= 6 x 10 ^ T^ 2 ergs/cm - sec - °K , (8) 


= 1.2 x 10 ^ T^/ 2 gm/cm 


sec 


(9) 


If fact, we are concerned here with a somewhat more complex hydrogen -helium 
plasma, and it is shown in Appendix I that the viscous coefficient is signi- 
ficantly composition dependent. For a 10$ helium admixture Yl should be 
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reduced to 0 .69 % , and we shall refer to this value as the coronal coeffi- 
cient of viscosity; the thermal conduction coefficient is essentially inde- 
pendent of mixing ratio and "K can be set equal to 

It is generally assumed that the coronal flow appears to an Inertial 
observer to be stationary, radial, spherically symmetric, and unaffected by 
electromagnetic forces, and these simplifications are tentatively adopted 
here; we also assume that no external heat source is present in the corona 
so that ^g is simply the acceleration of solar gravity, [as in (i), Eqs. 

(l) — (3) are to be integrated numerically and any discrepancy between 
the observed and predicted corona and solar wind will be' taken to indicate 
that at least one of the above assumptions is invalid. J For this case, 

the momentum and energy equations are, respectively. 


Nmu 


du 

dr 


+ £T (NkT) + 


NGM m 
® 

2. 


= it r &L ± (U) + iLry A 
3 dr dr V 3 Vdr 


1 d / 2 v 
(r u)l 


and 


SSL. . ?JSL. + 1 kT . X(T)r 2 dT _ k (T)r 2 ( du 
2 r 2 c dr 3 c l u dr 



( 10 ) 

( 11 ) 


2 

where E and c = Nur are constants of motion. It is convenient to introduce 
dimensionless parameters ^ = mu /kT Q , A = GM^m/kT^r, t = t/T q , e = E/kT Q 
(T is an arbitrary reference ' temperature ) , and Eqs. (10), and (ll) may then 
be rewritten in the form 


1 

2 


(1 - 


1 'i i 2x 

V 1 dA " 1 " A " 


dr 

5K 




1 1 /di\ 2 

WVa ) ‘ 


dA 1 


2iy 

A 2 



(12) 
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l-M^' 5/2 iK-< 5/2 (iKi)- • <«> 


Here 


A = 


2K(T o )GM c m 

k 2 T c 
o 


2 W k 

' B 'l Wi A = 0 - 69B c 


3. THE LOWER CORONA 

Since B is very small, (B/A = .0123 for a 10$ helium corona) the viscous 
terms should be negligible wherever the velocity gradients are "moderate", and 
Eqs. (12), (13) become 


1 /, _t_ \ d4^_ ^ , 2 t dr 

2 u • H> > dA '• 1 ■ 7\ ' 7i 


2 



(14) 

(15) 


In (i) it was assumed that this approximation is valid everywhere and the 
equations were integrated on an IBM 7090 computer inward from the earth to 
the sun with N e (l AU), u(l AU) and T(l AU) fixed at 3-4 cm”^, 3$2 km/sec and 
2.77 x 10 ^ °K, respectively (A = ^00, T q = 2 x 10^ ok). The thermal gradient 
at the earth was varied until a solution which has as >A C was ob- 

tained. A power series expansion about this singularity was then used to con 
tinue the integration into the subsonic lower corona |in (i), r Q = 5 >25 R @ , 

T c - H' c * 0.586]. 


The. density distribution predicted in this manner is the solid curve 
labeled "solar wind solution" in Fig. 1, and the dashed curves represent 
observed equatorial densities (the complete references to the experimental 
curves are contained in I. In Fig. 3 of Ref. (i), the measured electron 
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densities were compared with the computed total particle density, N = 1.9 N e 
for 10$ helium. The figures in the present paper contain more meaningful 
comparisons since all densities refer to electrons only. ) 

The preliminary solar wind density curve in Fig. 1 obviously differs 

from the experimental curves by at least a constant factor, but the general 

agreement in slope is reasonable. This must be contrasted with the behavior 

of the curve labeled "subsonic solution". Here N , T and dT/dr at r = 1.25 R 

are fixed at values appropriate for the solar wind expansion, but u(l.25R ) 

© 

is reduced from its critical or solar wind value I u (1.25R ) = 9*1^ km/sec I 

* s *w © ^ 

to 9 .04 km/sec. As the fluid moves outward, the new velocity distribution 
then remains just below the wind profile out to the = T c crossover at 
r = 5*25 R ; however, in the subsonic case u(r) then falls very rapidly 

C 0 

O 

as r is increased and since Nur is constant, the subsonic N(r) curve 
develops a very sharp knee near r = r Q . Thus, the smooth electron density 
profiles measured in the lower corona (say, r ^ 10 R @ ) already suggest very 
strongly that the corona flows continuously in the solar wind mode, rather 
than in any of the subsonic patterns. [_The term "subsonic" is used to 
designate solutions for which ^<t in the entire lower corona, but we do 
not exclude the possibility that t becomes equal to r at sane great dis- 
tance from the sun. The critical solution of Fig. 1 has e = 7*826, and 

s * w • 

-5 

the "subsonic" solution has e=€ -2.8x10 , thus the latter is not a 

s.w. 

solar "breeze" solution with e = 0. However, as /\— >0 (r — »oo ), Eq. (15) 
yields ^ /2 + (A/7) d(r^ 2 )/d)\ - e — > 0 so that, we might conceivably have 
*V(?\ = 0) = 0, t — * c 7\ 2 / T with Ac^ 2 = Je, in which case the solution 
would indeed be subsonic everywhere. On the other hand, a very small change 
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and in 


in the initial conditions then yields ^ —fc'V , 7( e _ 1^/2) = Ac"^ 2 , 
this more general case ^ becomes greater than t as r— *oo; at great dis- 
tances the apparently subsonic solution to Eqs. (l4), (15) with 0 < e <. e g 
generally develops into a supersonic wind.J 

The computed density profile in Fig. 1 is consistently below the ob- 
served distributions, and it is of interest to vary the theoretical parameters 
in order to determine the best fit. The parameters which might be changed are 
the conduction coefficient A, the mean mass m, the reference temperature T q , 
and the coordinates of the crossover point, J\ c and . 

The relative insensitivity of the predicted density profile to variations 
in composition (m) and base temperature was discussed in Ref. (i), and it was 
noted that the observed densities within 2 R^ were considerably higher than 
any computed ones. Table 1 shows that N g (r) is also roughly independent of 
the crossover location. Here N 0 (r) is computed for a range of ^ ~ h c 
combinations with A, T q and m fixed at 200, 1.5 x 10^ °K and 0.62 m^, res- 
pectively. The important result is that although N g does vary somewhat as 

^ and }\ are changed, the observed electron densities are always consider- 
c c 

ably greater than the predicted values, especially in the region below 2 R • 

© 

The main variation in the computed density profile is related to a change 
in (Nur ), or alternatively a change in the conductive coefficient, A. In 
Fig. 2 the extent to which the N g (r ) - curves depend on this parameter is 
shown There = O.586, and = I.365 (A = 100, 400, 600), )\ = I.32 
(A = 200); T q has been adjusted to give the best fit in each case, and all 
curves correspond to a uniform 10$ helium concentration! . For A = 100, 
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T q = 1.5 x 10 6 °K the conductive heating calculation yields an N g (r) distri- 
bution which agrees extremely well with the values obtained by Blackwell 

(1956) and Miehard (195*0 in the region 2R < r < (8 - 10) R . Beyond 

© © 

(8 - 10) R^ the predicted distribution is somewhat lower than that given by 
Blackwell, but it is known (Blackwell and Ingham, 1961) that the interpre- 
tation of these observations is unreliable beyond about 8 R @ ; it has been 
suggested (Chapman, 1961) that the densities beyond 8 R should be reduced, 
with N e (20 R @ ) shifted to perhaps one quarter of the value shown in Figs. 

1, 2. With these comments in mind, it may be stated that there is no de- 
finite discrepancy between the distribution predicted by the conductive 
heating theory and observations in the range 2 R @ < r < 20 R @ . However, 
if the reduction factor of four is meaningful, it seems likely that 
Blackwell's data should be reinterpreted so that N s (r) is less than 

S """" r 

N S * W * (r) (A = 100, T = 1.5 x 10 6 °K) for 8 R < r < 20 R . This would 
indicate that the solar wind actually flows more rapidly than predicted by 
the conductive heating theory beyond 8 R q . 

More positive conclusions can be reached concerning the discrepancy 

which is found in the shell between the coronal base (r — 1.06 R ) and 

© 

(2 - 2.5) R @ . Here the density observations are quite accurate, and all 
parameters (A, m, T q , and l^) have been varied in an attempt to re- 
produce the observations. However, the predictions based on Eqs. (1*0, 

(15) cannot be made to agree with the experimentally determined values, 
and we conclude that Eqs. (10), (ll) (viscosity is unimportant here) are 
not valid below (2 - 2.5) R^. The simplest possible explanation for this 
breakdown is that the external hydromagnetic or mechanical energy supply 
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which maintains the coronal base at a temperature near 1.5 x 10^ °K is 

actually dissipated over the large shell extending between 1.0 6 R @ and 

(2 - 2.5)R . In this case, Eqs. (10), (ll) are incorrect because the 
© 

presence of an important local power source is omitted. However, while 
it is undoubtedly true that the external heat source penetrates some 
finite distance into the corona, several other reasonable processes can 
invalidate Eqs. (10), (ll) in this region, and at present it is impossible 
to identify the thickness of the heating shell with confidence. For in- 
stance, it has been suggested (L. Davis, Jr., private communication) that 
the filamentary structure and other deviations from spherical symmetry 
(sunspot associated arches, streamers, etc.) can have an extremely im- 
portant effect on the flow patterns near the coronal base; outgoing 
plasma streamers may be accelerated as they flow through magnetic nozzles 
associated with filaments, and it is therefore unlikely that the base 

region (r < 2 R ) can be well described hy Eqs. (10), (ll), even in the 
© 

absence of an external heating source. The flow pattern may also be 
distorted by diffusive settling of the helium component leading to a 
variation in the mean mass with altitude (Parker, 1963). Another cor- 
rection which might be significant is associated with the weakening of 


the solar gravity field by centrifugal forces. Near the base of the 

corona the flow is not radial and uv(r-* R )~K£LR sin 0, with 

t> © © ’ 

2 . 94 x 10 rad/sec. It is not known just how \x^ decreases with 

increasing r, but even during quiet periods ^(r > R^) is certainly 

greater than the angular momentum conservation limit, Hr? sin 9/r, be- 

© 

cause of the interaction between the plasma and the general magnetic field 
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(Parker, 1958; Chapman, 1961; Scarf, 1964); thus, the effective solar gravity 
field can be distorted by corotation phenomena over a range of several solar 
radii. On the basis of analysis of Pioneer 5 field and Mariner 2 plasma 
measurements, Greenstadt (1964) found some observational support for a 
strong corotation model. The second pair of Vela satellites carry plasm 
probes with high angular resolution and the preliminary solar wind obser- 
vations also suggest a significant amount of corotation (S. J. Bame, Pri- 
vate communication). It is likely that all of these effects contribute 
to the discrepancy near the coronal base. 

4. VISCOUS EFFECTS 

The search for a best with the conductive heating equations 
involved adjustment of parameters and numerical integration in the region 
R g) <r<20 R^. More explicitly, boundary conditions were fixed near the 
sonic transition fr Oi (5 - 7 ) R for the curves in Fig. 2 J , and the 
numerical integration of Eqs. (l4), ( 15 ) was carried out in both direc- 
tions away from the crossover. It might be assumed that as tends to 
zero, t then automatically vanishes and ^ tends to a constant asymptotic 
streaming energy. 

In fact, the numerical integrations do not always yield the temperature 
and kinetic energy profiles which are anticipated on simple physical grounds. 
Typical numerical results are shown in Fig. 3> and it can be seen that 
Tg-40 for 7\ >0, and T^(o) > 0. These anomalous examples have straight- 
forward analytical expressions; for small A the conductive heating equa- 
tions, (l 4 ), ( 15 ) possess supersonic solutions with a variety of asymptotic 
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forms, some of which are physically unreasonable. 


By inspection of Fig. 3 > it can be determined that the solution with 

vanishing temperature at finite A = A has t,— *( A~A ) 2/ ^ 7 . This unphysi- 

O Ct o 

cal behavior actually does represent a true numerical solution of Eqs. (l4), 
(15); since the expansion 

t-4c() - TO\-\) 4/7 /3Ac 3/2 + ••• 
ty-*2(e - Ac 7 ^ 2 /7) + 2 A - 2 t + ••• 

formally satisfies the coupled equations for all c. Similarly, the curves 
with finite temperature at infinity appear to correspond to the set with 
asymptotic solutions 


[ A-»A 0 > 


r 7/Z -> * \ fj? A ♦ d 1 * »„/*> > in > + 

V— >c(e,T a> ) - 4 T^.fi.n A + ••• 


A->o. (17) 


However, the continuous family of physical solutions must have t( 7\= 0) = 0, 
and this generally involves 


c> 2 / T + (77/9Ac 3 ^ 2 ) A 4 ^ 7 + •• 
^-»2(e - Ac 7//2 /7) + 2 A - 16 t + •• 


A->° 


(18) 


We regard Eq. (l8) as the acceptable physical family of solutions to 
describe an inviscid fluid with streaming and conductive energy transport 
which both remain finite at infinity. Parker's analytical studies of Eqs. 
(l4), (15) (1964 a, b; 1965) have been based on the assumption that Eq. (l8) 
with arbitrary c is the correct description of the flow in the small A 
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region. In fact, Eqs. (l4), (15) admit another special asymptotic behavior 
for the particular case in which the asymptotic energy transport is purely 
kinetic (Whang and Chang, 1964); i.e.. 


Y- 


(35>/2A) 2/5 + ... 


’ 2e + 2 > - 12t + 


A->o 


(19) 


This overabundance of asymptotic forms indicates that difficulties can 
be encountered when a numerical integration of Eqs. (l4), (15) proceeds in 
the direction of decreasing 7i • However, Eqs. (l6) - (19) illustrate a 
more significant physical fact; since Eqs. (l4), (15) possess unreasonable 
solutions, the equations themselves cannot represent an adequate description 
of a physical fluid. 


It is not difficult to identify the unphysical aspect of Eqs. (l4) and 
(15) because no real fluid can have a vanishing Prandtl number. The exis- 
tence of anomalous solutions to the inviscid equations indicates that 
viscous effects must indeed be significant so that the full Navier-Stokes 
equations (12, 13) should be retained. 


By attempting to treat the viscous terms in Eqs. (12, 13) as a small 

perturbation, we have determined that these corrections do remain negligible 

below the sonic transition point, (i.e., A>A )• At the viscous 

c c 

contribution to the energy equation begins to be significant (on the order 
of a few percent of the thermal energy), and we therefore carry out all sub- 
sequent calculations using Eqs. (l4), (15) for h>A\ -S ) and Eqs. (12), 
(13) for J\< ( A c - S )j the results are not sensitive to the precise value 
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for ^ if <fc/K>« 1, and this approximation allows one to use the 
simple techniques developed in (i) for treating inviscid flow near the 
sonic transition. 

The dimensionless viscous coefficient appropriate for a 10$ helium 
corona with A = 200 is B = 2.46. The heavy curve in Fig. 4 shows the 
streaming energy profile obtained by numerical integration of the complete 
hydrodynamic equations, including uninhibited viscosity beyond the transi- 
tion near 7*5 R (T * 1.5 x 10^ °K). The dashed curve with B = 0 is the 
© o 

original solution without viscosity; and the temperature curves are es- 
sentially independent of B in this region. The most striking aspects of 
these calculations are (l) the fact that the streaming energy profile is 
greatly changed in magnitude by the inclusion of viscosity with 

u^(20 R , B = 2.46) > 2 u^(20 R , B * 0); (2) the observation that the 
© © 

viscous term appears to operate in the wrong direction by yielding a local 
speed which is greater, not less, than that obtained when viscosity is 
ignored. 

In order to investigate the origin of the second effect, numerical 
integrations of the Navier -Stokes equations were repeated using artificial 
enhanced and inhibited values for the viscous coefficient (B = 0.1, 1.0, 
3.56). The effect of this variation in B is also shown in Fig. 4 (all 
other parameters and boundary values remain the same). It can be seen 
that as B is increased from 0.1 to 3*56* the streaming energy at any given 
radius is reduced as might be anticipated. Larger viscous dissipation 
does indeed produce lower flow speeds relative to the sun. However, it is 
also clear from Fig. 4 that the flow pattern with B = 0 is completely 
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unrelated to that obtained in the limit B — - ^ 0 . 

From a mathematical point of view, the fact that B s 0 is a singular 
limit is not difficult to understand. The coupled Navier -Stokes equations 
are solved by formally writing Eq. (13), for instance, in the form 



U dr/d , A, A, e) 

B 


(13a) 


and this display already casts doubt on the relevance of the inviscid energy 
equation. A similar situation occurs when one examines the limit A— >0 for 
the inviscid model. For instance, Eq. (15) has the form 


dT T,>, €) 

d>> . 


and the acceptable solutions j^Eqs. (18), ( 19 )J have temperature expansions 
in A” 1 which obviously become meaningless in the limit A—>0; these con- 
ductive heating solutions are not related to the bounded adiabatic solution 
of Eqs. (lk), (15) with A = 0. By the same token, it appears that Fig. k 
merely illustrates that fluid flow with a small amount of viscous dissi- 
pation differs greatly from inviscid flow. This conclusion is a familiar 
one to workers in aerodynamics. 

We adopt the solutions with B/A = .0123 as "the correct ones for the 
solar corona in the regions where the Navier -Stokes equations are applicable. 
The related increase in the local streaming speed requires a corresponding 
reduction in N with N(r,B) u(r,B) = N(r,0) u(r,0) since the same value for 
the constant (Nur ) was used throughout. This means that the predicted 
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values for N (r) should he lower than those shown in Fig. 2. At r — (7-10)R 
6 0 

no significant reduction is needed hut N e (r,B)/N e (r,0) — 0.77 at r = 15 R @ , 
0.6l at 20 R^, etc. These changes increase the discrepancy between the pre- 
dicted Ng(r) curves and Blackwell's original results in the uncertain region 
(10-20 )R @ , however, the best theoretical curve is now brought into closer 
agreement with the corrected density profile discussed by Chapman. It is 

hoped that better density measurements for the region (10-20 )R will be 

0 

available in the near future so that more meaningful comparisons can be made. 
5. CONCLUSIONS 

It has been shown that in the region between 2 R and 15 R , where 

© © 7 

reasonably accurate density measurements are available, the solar corona 

can be well described by the conventional Navier -Stokes equations with no 

external heat source. The best fit to N 0 (r) is obtained for a 10$ helium 

corona with T q — 1.5 x 10 °K, A = 100, B = .0123* Some representative 

numerical values for this solution are given in Table 2, and it can be seen 

that the streaming speed increases very rapidly beyond the sonic transition 

near 7 R • This strong acceleration in the region (10-15) R is a consis- 
0 0 

tent feature of the full Navier -Stokes equations, but it is not obtained 

for the singular and unphysical case of zero viscosity. Beyond 15 R , the 

© 

velocity gradients are so steep that the Navier -Stokes formalism becomes 
questionable. This problem will be discussed in more detail in the final 
paper of this series. 

Below 2 R there appears to be a genuine and significant discrepancy 
© 

between the predictions of the idealized theoretical models and observation. 
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This suggests that one or more of our assumptions (spherical symmetry, uni- 
form composition, no corotation, no external heat source) is invalid in the 
coronal base region. 
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APPENDIX A 


Viscous Coefficient for a Hydrogen -Helium Mixture 


The viscous transport coefficient for fully ionized hydrogen has been 
evaluated by Chapman (195*0, who points out that the electrons do not 
appreciably affect the viscous energy transport due to their small mass. 

The appropriate moments of the Boltzmann equation for a mixture of fluids 
are discussed by Burgers (i960). In this appendix Burgers' equations are 
applied to a mixture of 90 $ hydrogen - 10 $ helium. The effect of electrons 
is not explicitly included. 


Burgers' Eq. (5 - 3*0 is Just the i-j moment of the Boltzmann equation 
in concenient form. The usual linear relation between viscous stress and 
strain is obtained by neglecting time derivatives and gradients; the re- 
sulting equations for a proton -alpha particle mixture are the following: 
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Here the subscripts p and a refer to protons and alpha particles respectively, 
P^ and are the viscous stress and strain tensors, respectively, and 
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The Z functions are appropriate integrals of the collision cross sections; 
they are evaluated by Burgers for Coulomb interactions. The following 
particular integrals are useful here: 

Z pp l} 2Z =2Vn , (^7f)logA ; A = T § — 1 
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Inserting these expressions, one can solve Eqs 0 (A-l) and (A-2) simultaneously. 
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where f = — , 
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The usual coefficient of viscosity for such a mixture is defined by 
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One can rewrite Eq. (A-l) as follows: 
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For a pure hydrogen plasma one has simply 
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Combining Eqs. (A-k), (A-5) and (A-6) yields finally 
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For the case of interest here, in = km and f = From Eq. (A-3) one 

Cc p y 

then obtains A(f) = O.O89. Inserting these values in Eq. (A- 7 ) yields for a 
90% hydrogen - 10$ helium mixture the result = 0.6$Y^ o * 
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TABLE 1. The variation of the density profiles with changes in 

2 

crossover location (A c = GM^m/kT Q r c ) and energy (v^ c = mu c /kT Q ) . 

_ O 

The tabular entries represent electron density (cm J ) and these 

results are all derived using Eqs. (l4), (15) with A = 200, 

m = 0.62 m and T = 1.5 x 10^ °K. 

P ° 



r= - 

0.585 

A c ■ 1-365 

r/R 

A c = 1.32 

\ = 1.40 

,x c 

^> c = 0.54 

V|/ c = 0.62 

1,6 

4.61 x 10 6 

3.38 X 10 6 

3.49 X 10 6 

4.40 x 10 6 

2.0 

1.39 x 10 6 

1.09 X 10 6 

1.16 x 10 6 

1.32 x 10 6 

4,5 

5.45 x 10^ 

4 

5.05 x 10 

4 

5.31 x 10 

4 

5.20 x 10 

6,8 

1.58 x 10^ 

1.51 x 10 4 

li 

1.60 x 10 

4 

1.50 x 10 

7 2 

1.34 x 10^ 

1.28 x 10^ 

4 

1.36 x 10 

1.27 x 10 4 
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TABLE 2. The best fit solution to the Navier -Stokes equations in 

the region 2R < r < 15R . This solution (m = .62 m , A = 100, 

© © P 

6 

B = 1.23, T = 1.5 x 10° ok) has a sonic crossover near 7*25 R • 
o © 

Viscous effects are significant only beyond the crossover radius 
and the Navier -Stokes equations become questionable beyond 15 R . 
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N (cm" 3 ) 
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18.4 

2.03 x 10 6 

1.47 x 10 6 

4.5 

72.4 

1.02 x 10 5 

1.03 x 10 6 
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1 — 1 

3.09 x lO^ 

8.89 x 10 5 
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no 

2.62 x 10 k 

8.73 x 10 5 

10 

146 

1.02 x 10 1 * 

7.92 x 10 5 

12.6 

190 

4.95 x 10 3 

7.55 x 10 5 

15 

246 

2.70 x 10 3 

7-35 x 10 5 
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FIGURE CAPTIONS 


Figure 1. Comparison between sample predictions of the conductive heating 

equations (no viscosity) and observed electron densities in the 

inner corona. The references to the experimental papers are 

contained in (i). Hie theoretical curve labeled "solar wind 

solution" represents the density profile for flow which becomes 

supersonic in the neighborhood of 5R . The "subsonic solution" 

© 

then represents flow in which the energy per particle is reduced 

5 

by about one part in 10 ; beyond 5R the subsonic speed falls 

© 

rapidly so that the density becomes large. 

Figure 2. Effect of varying A and T q in the conductive heating equations. 

Beyond about 8R @ the A = 100 predicted density is somewhat below 
the observed curve (Blackwell, 1956) but the interpretation of 
the measurements is not reliable in this region. The discre- 
pancy for r< (2-2.5)^ probably represents the appearance of 
complexities which invalidate the model, as described in the 
text. 

Figure 3 * Typical anomalous solutions of the inviscid model equations. 
These solutions are described by Eqs. (l6), (17). 

Figure 4. Effect of viscosity on coronal flow. The heavy curve with 

B = 2.1^6 (A = 200) represents the streaming for a 10$ helium 
corona. The dashed curve (B = 0) is the corresponding solution 
for a hypothetical non-viscous fluid; the temperature is essen- 
tially independent of Y^(or B) in this region. The solutions 
for artificial values of viscosity (B = 0.1, 1.0, 3.56) are 
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included to illustrate that the limiting flow for small B is 
unrelated to the singular result obtained with no 


viscosity. 
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